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FOREWORD 


This report, prepared by the Martin Marietta Corporation, 
Denver Division, under Contract NAS8-30690, presents the results 
of an analytical and experimental study of transient liquid 
motion similar to that encountered in orbiting spacecraft. The 
study was performed from March 1974 to February 1976 and was 
administered by the National Aeronautics and Space Administration, 
George C. Marshall Space Flight Center, Huntsville, Alabama, under 
the direction of Mr. Frank Bugg. 

Some of the results of this study were published in an 
interim report, MCR-75-4, February 1975. In addition, two 16 mm 
high speed movies have been generated which document the test 
results . 
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ABSTRACT 


This report presents the results of a two year study of 
transient liquid motion such as that which will be experienced 
during orbital maneuvers by Space Tug. A test program was con- 
ducted in the Martin Marietta Drop Tower Test Facility involving 
forty-five drops. Biaxial, low-g accelerations were applied to 
instrumented, model propellant tanks during free-fall testing, 
and forces exerted during liquid reorientation were measured and 
recorded. High speed photographic records of the liquid reorien- 
tation were also made. The test data was used to verify a 
mechanical analog which portrays the liquid as a point mass 
moving on an ellipsoidal constraint surface. The mechanical 
analog was coded into two Fortran IV digital computer programs: 
LAMPS; Large Amp litude ^losh (a two dimensional simulation), and 
LAMPS3 (a three dimensional simulation, cast in the framework of 
a general spacecraft simulation program) . Results show excellent 
correlation between test data and analytical predictions of 
reorientation forces and liquid center of mass motion, verifying 
the basic analytical approach. 
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I. INTRODUCTION 


In the near future, vehicles similar to Space Tug will 
perform orbital maneuvers while carying a large mass of liquid 
propellant. An indepth understanding of the interaction forces 
between the propellant and the space vehicle is required to 
properly assess the dynamics of these maneuvers. During orbital 
maneuvers, the propellant mass is subjected to small accelera- 
tions which can induce large amplitude propellant slosh. Due to 
the relatively large mass of propellant, the forces exerted on 
the spacecraft by the moving propellant may have a significant 
effect on gross vehicle motion. Knowledge of these interaction 
forces is imperative in the design of spacecraft control systems, 
propellant management devices and docking mechanisms, in addition 
to being required for maneuver studies such as deorbit/reentry, 

A two year study (Phase I - 1974, Phase II - 1975) has been 
conducted to develop and verify a mechanical analog to simulate 
large amplitude liquid motion in a container subjected to low-g 
acceleration environments. The primary purpose of the model is 
to simulate the interaction forces between the liquid and its 
container. The study consisted of both experimental and analyti- 
cal tasks. 

The Martin Marietta (Denver) Drop Tower Test Facility was 
used in the conduct of the experimental phase of the study, A 
test module capable of measuring forces generated by large ampli- 
tude slosh in a scale model tank was constructed. Dimensional 
analysis techniques were used to ensure that the model simula- 
tions would be representative of reorientation in a full scale 
liquid oxygen tank. Three model tanks were tested; all three were 
axis 3 mimetric with cylindrical barrel sections and hemispherical 
domes. The models consisted of a baseline tank, a half size tank 
(to increase equivalent reorientation time), and a baseline tank 
containing three ring baffles. During the study, the test module 
was dropped in the free-fall tower (simulating low-g) and small 
biaxial accelerations were applied. The ensuing liquid motion 
was photographed and two dimensional forces were measured and 
recorded. A total of forty-five tests were conducted. Various 
tank fill volumes, tank orientations and acceleration magnitudes 
were investigated. The test durations (1.5 to 2.0 seconds) 
correspond to approximately 15 to 20 seconds of liquid motion in 
a full size tank (10 foot diameter). Chapter II details the 
test program. 
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A two dimensional mechanical analog was developed during 
Phase I to simulate the observed large amplitude slosh. The 
analog portrays the liquid as a point mass moving on an ellipsoi- 
dal constraint surface. The constraint surface is defined by the 
ellipsoid which best fits the locus of liquid center of mass 
locations prescribed by slowly rotating the tank in a one-g field. 
The mechanical analog was implemented in a computer program, LAMPS 
(Large Amp litude ^losh) , which outputs force time histories on the 
tank due to liquid motion. The Phase I test program and computer 
program LAMPS were both detailed in an interim report (Reference 
1 ). 


During Phase II of the study, the mechanical analog was 
expanded to three dimensions and integrated into the general 
spacecraft equations to facilitate its use in subsequent analyses. 
The three dimensional model was implemented into a computer pro- 
gram LAMPS3, which is in the framework of an existing dynamic 
simulation program DYNAMO (Dyn amic ^alysis Multiple ^tions) . 
Chapter III discusses the implementation of the model; Appendix A 
presents a users guide for LAMPS3. 

Chapter IV presents a discussion of the acquired test data 
and correlation with analytical results. 
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II. EXPERIMENTAL INVESTIGATION 


The primary objective of the experimental investigation was 
to generate data for correlation with the computer model developed 
under the analytical task. The experiments simulated the reorien- 
tation of liquid within a propellant tank as it would occur during 
spacecraft maneuvering. The motion of the propellant in the main 
tanks of a Space Tug during docking maneuvers was used as an 
example to define typical environmental conditions. Scaling was 
used to relate those conditions to a small tank model. 

A drop toxv^er, which can simulate a lox> 7 -g environment, was 
selected as the means of performing the tests. The testing was 
conducted in tx^^o phases, x^ith 22 tests in Phase I and 23 tests 
in Phase II. For each test the motion of the liquid was photo- 
graphed and the forces exerted by the liquid on the tank were 
measured. Since the test conditions used have not been experi- 
mentally simulated before, the tests add to the basic understand- 
ing of the reorientation of propellant within a tank. 


A. TEST APPROACH 

When the propellant within a tank is reoriented, it is 
displaced from its initial position and relocated to a new 
equilibrium position due to an acceleration field acting on the 
tank. A prior maneuver is assumed to have positioned the liquid 
at its initial location. Since an attitude control system is 
defined as the source of the accelerations, the orientation of 
the acceleration with respect to the initial liquid position is 
arbitrary. For the test program, the most significant liquid 
motion and forces were desired, so the reorientation produced 
by an acceleration essentially opposite in direction to the one 
that established the initial position of the liquid was con- 
sidered. An off-axis component was added to the reorienting 
acceleration to avoid symmetrical liquid motion, which tends 
to be a special case. 

The typical conditions for a test are as shown in Figure 
II-l. A cylindrical tank with hemispherical end domes and an 
L/d of 1.28 was tested. A volume of liquid was initially at 
rest and the liquid surface was positioned with respect to the 
tank axes as determined by the angle 0x. The axial (Aa) and 
lateral (Al) components of the reorienting acceleration are 
defined with respect to the initial liquid interface. The axial 
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Aa 


Initial Liquid 
Interface 


Rgure IN. Initial Conditions 
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component of the acceleration is always defined as being 
perpendicular to the liquid interface. The acceleration was 
continuously applied and the liquid was reoriented to the oppo- 
site end of the tank. 

A scaling analysis was performed to determine the appropriate 
test parameters based on a Space Tug tank; 3.1 m (10 ft.) diameter, 
maneuvering accelerations on the order of 0.02 g to 0.04 g., and 
liquid oxygen as the propellant. A complete discussion of the 
scaling approach can be found in the Interim Report (Reference 1) ; 
only a brief summary will be presented here. 

The following dimensionless parameters characterize the 
motion of the liquid as it is reoriented. 

Froude Number - Ratio of inertia to gravity force 

Bond Number - Ratio of gravity to surface tension force 

Reynolds Number - Ratio of inertia to viscous force 

Froude number has been expressed as a function of the Bond 
and Reynolds numbers based on the results of numerous reorienta- 
tion tests. If both the Reynolds and Bond numbers are suffi- 
ciently large, which was the case for the selected test conditions, 
the Froude number is a constant. By equating the Froude number 
of the full size system and the model, similarity of test condi- 
tions is obtained. 

This scaling analysis yielded a relation between time (t) , 
acceleration (A) , and tank radius (r) : 

ta _ / Am ^ a , 

tm " W Aa r 

1 m 

where the subscript **a” refers to the full-size system and the 
subscript *'m” refers to the model. 

The purpose of the scaling was to make the test conditions 
representative of the actual operational environment. The test 
conditions were duplicated in the analytical model so scaling 
did not enter into the correlation. 

A model tank, 12.7 cm (5.0 in) in diameter, was used for 
the Phase I tests. For Phase II, this tank was again used in 
addition to a tank of the same size with ring baffles, and also 
a 4.6 cm (2.5 in) diameter tank. The dimensions of the three 
model tanks are shown in Figure II-2. 




Figure 11-2. Model Tank Dimensions (all dimensions are in inches) 
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The drop tower was selected as the means of simulating the 
acceleration environment. Following release of the drop capsule, 
a near zero-g test period of up to 2.1 seconds is available. 

During the test period, an acceleration can be applied to the 
model to simulate the reorientation acceleration. Due to physi- 
cal constraints of the drop tower, the test acceleration cannot 
be larger than about 0.1 g. 

A large model acceleration and a small model radius are 
desirable, based on the Froude number scaling, to make test time 
represent as long an actual time period as possible. A maximum 
model acceleration of 0.1 g was selected and, to consider the 
effect of acceleration, a smaller value of 0.05 g was also used. 
The lateral component of the acceleration ranged from a momentary 
pulse to a value of 50% of the axial acceleration. 

Too small a model radius would reduce the magnitude of the 
liquid forces, hence two tank sizes were considered. The large 
tanks favored the liquid forces and the small tank provided a 
longer full-scale test time. 

FC-43, a very dense fluorocarbon solvent, was selected as 
the test liquid. Its properties at 20°C (68°F) are as follows 
(Reference 2) : 

Density 1.905 gm/cc (118.9 lbm/ft5 

Surface Tension 16.7 dynes/cm (1.14 x 10 Ibf/ft) 

Viscosity 6.5 cp (4.36 x 10 Ibm/ft-sec) 

This choice of test liquid, tank size and accelerations 
yielded a time scaling ratio (ta/tm) of 7.4 for the larger 
tank. Hence, 2.1 seconds of test time were equivalent to 15 
seconds of real time. For the small tank, the ratio was 10.4 
and 2.1 seconds of test time equals 21.8 seconds of real time. 

Earth gravity was used to establish the initial position of 
the liquid in the tank. A flat interface is a realistic initial 
condition for a large tank in a relatively low-g environment. 

The acceleration would have to be 10 g or less with a tank on 
the order of 3 meters (9.8 feet) in diameter before the inter- 
face would have any significant curvature. 

The parameters for each of the tests are listed in Table 
II-l. The tank orientation is the angle 0 x defined on Figure 

Spring motors were used to produce the acceleration and the 
motor force is listed in the table. Elsewhere in this report, 
the actual acceleration achieved is presented along with the 
force data (see Appendix B) . 
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TABLE II-l TEST MATRIX 



LIQUID 

TANK 

SPRING MOTOR FORCE 


TEST 

VOLUME 

ORIENTATION 


N (LBF) 



NO. 

(%) 

(6X DEGREES) 

AXIAL LATERAL 

TANK 



PHASE 

_I 




1 

25 

0 

67 

(15) 3.3 

(0.75) 

Large 

Tank 

2 

50 

0 





3 

75 

0 






4 

25 

30 






5 

50 

30 






6 

75 

30 






7 

25 

60 






8 

50 

60 






9 

75 

60 






10 

25 

90 






11 

50 

90 






12 

75 

90 






13 

25 

0 

133 (30) 




14 

50 

0 






15 

75 

0 






16 

25 

45 






17 

50 

45 






18 

75 

45 






19 

25 

90 






20 

50 

90 






21 

75 

90 





1 

22 

10 

0 








PHASE 

II 




1 

25 

0 

133 (30 0 


Large 

Tank 

2 

25 

45 


0 


Large 

Tank 

3 

10 

0 


3.3 

(0.75) 

Large 

Tank 

4 

10 

45 




with 

Baffles 

5 

10 

90 






6 

25 

0 






7 

25 

45 






8 

25 

90 






9 

50 

0 






10 

50 

45 






11 

50 

90 






12 

10 

45 

67 

(15) 




13 

25 

45 

y 





14 

50 

45 

1 





15 

25 

0 

133 

(30 


Small 

Tank 

16 

25 

45 






17 

50 

0 






18 

50 

45 




f 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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TABLE II- 1 (Cont) 



LIQUID 

TANK 

SPRING MOTOR FORCE 


TEST 

VOLUME 

ORIENTATION 

N (LBF) 



NO. 

(%) 

(0X DEGREES) 

AXIAL LATERAL 

TANK 



PHASE II (Cont) 



19 

25 

45 

67 (15) 3.3 

(0.75) 

Small 

Tank 

20 

50 

45 

67 (15) 3.3 

(0.75) 



21 

25 

45 

133 (30) 6.7 

(1.5) 



22 

50 

45 

1 6.7 

(1.5) 

1 

f 

23 

25 

90 

T 3.3 

(0.75) 

f 
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B, TEST SYSTEM DESCRIPTION 

A test system that can produce the required subscale model 
test conditions and measure the liquid forces was designed and 
built for Martin Mar ietta * sDrop Tower Test Facility. Flexibility 
to duplicate the varied test conditions, and sensitivity to record 
the small liquid forces, were the key requirements in designing 
the test system. 

1. Test Module - The test module consists of the tank, 
force measuring links and slider mechanism. This module is shown 
mounted on the drop capsule in Figure II-3. Figure II-4 and 
Figure II-5 present front and back views of the box in which the 
force links and tank are mounted. 

The model tanks were made of clear plastic. The flange 
around the tanks provided structural strength and permitted the 
tanks to be mounted at the proper angle within the force link 
yoke. 

Three force links, two vertical and one lateral, allowed all 
forces acting on the tank to be measured. Bearings at each end 
of the links permitted only forces along the link axis to be 
measured. The bearings that are mounted on the box are self- 
aligning. 

Three flexures, perpendicular to the plane of the force 
links, prevented any motion of the tank out of that plane. The 
spring constant of these flexures is small in comparison to the 
spring constant of the load cells and force links. Therefore, 
the effect of these flexures on the force sensed by the load 
cells was insignificant. 

The platform of the slider was attached to the rails with 
three linear bearings; one under the camera and two under the 
tank. A constant force spring motor provided the lateral accel- 
eration of the slider. An electric solenoid was used to release 
the slider at the beginning of the test. 

2. Drop Test Facility - The complete drop capsule is shown 
in Figure II-6. Due to the rather high accelerations being used, 
evacuation of the drag shield was not necessary. A simple frame 
was mounted over the test module rather than sealing the drop 
capsule with its cylindrical cover. The spring motors that pro- 
vide the axial acceleration of the drop capsule and a crush tube 
are mounted on the conical base. 


Lateral Spring Motor 
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Test Module Mounted on Drop Capsule 





Figure II-4. Camera View of Tank and Force Links 




Figure II-5. Back View of Tank and Fore 
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Complete Drop Capsule 
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The total drop test system is illustrated in Figure II-7. 
The cable from the axial spring motor was extended and secured 
to the bottom of the drag shield. After releasing the drag 
shield from the top of the 23-meter (75-foot) drop tower, a 
solenoid was actuated releasing the slider. Both the lateral 
and axial spring motors accelerated the test module throughout 
the drop test. The drop capsule impacts the drag shield, with 
the crush tube absorbing the impact, and the drag shield lands 
in a bin of wheat at the end. of the test. 

3. Instrumentation - The motion of the liquid was recorded 
with a 16-mm Milliken DBM-3a camera mounted on the slider. The 
film speed was 200 frames per second. Immediately before the 
drag shield was released, the camera was started and it was 
automatically stopped when the drag shield impacted the wheat. 

Quartz crystal load cells (Kistler Model 912) were used to 
measure the liquid forces. These load cells have a capacity 
of 2220N (500 Ibf) in tension and 22200N (5000 Ibf) in compres- 
sion, providing the capability of withstanding the impact at 
the end of the test. Peak, high frequency accelerations of up 
to 160g have been measured at impact. Due to their high degree 
of linearity, these load cells are fully capable of measuring 
the small forces due to the liquid motion. 

The load cells were mounted in the force measuring links. 
Low noise cables were used to feed the output of the load cells 
to charge amplifiers. The charge amplifiers were located about 
half way up to the drop tower to minimize the motion of the 
cable as the drag shield falls. The amplifiers were set on 
long time constant and the most sensitive scale that could be 
accommodated, to measure the low amplitude and low frequency 
forces. Each charge amplifier input was momentarily grounded 
prior to the test, so all forces were measured with respect to 
zero at one-g. 

The output of the charge amplifiers was fed in parallel to 
both a tape recorder and a chart recorder. In order to filter 
out the vibration induced by the camera motor, a 10 Hz low-pass 
filter was used in the amplifier for the chart recorder. An 
end-to-end calibration of the force measuring system was accom- 
plished with the fixture shown in Figure II-8. Known weights 
were suspended from the hook at various positions with respect 
to the force links, and the charge amplifiers were adjusted to 
give the proper output. 
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An accelerometer was mounted on the slider for the Phase I 
tests, to measure the axial acceleration of the drop capsule. 
It was found that the acceleration can be readily calculated 
from the known capsule travel distance, weights and test time. 
The accelerometer was not used during Phase II. The lateral 
acceleration can be calculated from the travel versus time 
obtained from the film data. 
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Instrumentation Cable 


Model Tank 


Slider 


Drop Capsule 


Axial Spring Motor 


Drag Shield 


Figure II-7. Conqjlete Drop Test System 
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Figure H-8. Calibration Fixture 
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C. DATA REDUCTION 

A lateral force (F ) and two vertical forces (F and F ) were 
measured and recorded during the drop tests (see Figure II-5) . 

The test data was manually scaled from oscillograph records and 
converted to a punched card data bank. A data reduction computer 
program was then used to convert the raw data from voltage to 
engineering force units using appropriate scale factors based 
on charge amplifier sensitivities. The data was then smoothed 
to remove test fixture and camera noise by a moving average 
digital low pass filter set at 10 Hz, Figure II-9 depicts the 
shape of the filter used. The force triad was then transposed 
to the tank triad as shown in Figure II-IO. The following set 
of equations was used to perform the transposition. 
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where subscript (I) denotes the intertial triad and subscript 
(T) denotes the tank triad. The results were plotted with time 
as the ordinate. 

To facilitate comparison between the test and analytical 
results, the test data was further adjusted. As previously 
mentioned, the force gages registered "0" in 1 g prior to each 
drop. The analytical model records this one-g force as a nega- 
tive force in the direction. To make the analysis and test 
results compatible, the initial zero test forces were converted 
to negative forces. In this conversion, the fluid cm offset 
was accounted for in the determination of the adjusted F„ and 
F values. This conversion allows direct comparison between 
predicted and measured force time histories. 





III-l 


III. MECHANICAL ANALOG 


This chapter details the approach whereby the tank/liquid 
mechanical analog is cast into a state space framework that is 
amenable to accommodation into a typical spacecraft system of 
governing equations. The basic approach is characterized as a 
system of interconnected "bodies" subject to certain constraint 
conditions which restrict the relative motion between the 
"bodies", one of which represents a typical liquid mass within a 
container. The ensuing discussion describes some of the particu 
lars relating to the tank/liquid mechanical analog which depicts 
the fluid mass as moving on a predetermined constraint surface 
within the tank boundry. An in-depth theoretical development 
for the general class of spinning, interconnected "bodies" can 
be found in Reference 3. This discussion will be limited to a 
cursory overview of the developments contained in the reference 
and focus on some of the specifics related to the three dimen- 
sional tank/liquid mechanical analog developed during Phase II. 


A. GOVERNING EQUATIONS - GENERAL 

The system of equations governing the dynamics of an 
mterconnected set of bodies undergoing large relative motions 
IS depicted in a state space form that can be readily solved by 
numerical mtegration^techniques on a digital computer. In terms 
of the state space, y , a canonical, first-order expression of 
the governing equations appears in a rather general form as; 


d i ^ / i . 

dt y = f (y . t) (iii-i) 

More specifically, the matrix form for the ‘equations as applied 
to our class of problems i s given as 

I I +[b]^ ) X 1 

{ ^ t ^ [b] j ( b (1II-2) 

[b])„| . 1»| 


c 
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State variables of the configuration space, y^, include 
ordinary momenta, p , and relative position and/or attitude 
coordinates, /S The vector, ^ contains such items as Euler 
angles and inertial position coordinates. The vector, \ j , 
represents the constraint forces and/or moments required to 
affect the particular constraints for the problem (e.g., hinge 
constraints between adjacent bodies and/or motion restricted to 
a constraint surface) . 


For a typical body, k, of the system, the component ordinary 
momenta vector, | P is 


f = ["]k i“i 


(III-3) 
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The matrix, |^t) J , contains kinematical coefficients which are 
described in further detail in the following subsection. The 
matrix, [o]. (for the kth body) contains three (3 by 3) parti- 
tions of skew symmetric submatrices of the body velocity vectors, 
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The matrix, LbJ. , contains kinematical coefficients of the same 
form as matrix, 3 . The vector, in the last expression 

of equation III-2 is referred to as a "rhenomic" prescription of 
relative velocities within the system. Whenever {a 0 , 

an ordinary constraint is imposed that fixes certain degrees of 
freedom with respect to each other within the overall dynamical 
system. The rhenomic constraint (time dependent) results from 
prescribing the rhenomic velocities , ^ ^ , as a function of the 

independent variable, time. This latter situation will be 
employed to impose the constraint whereby the tank fluid mass is 
restricted to move on an ellipsoidal surface (taken with respect 


to the tank axis system). There exists 
formation that relates the nonholonomic 
generalized velocities , 


another important 
velocities, fu/ , 

K. 


trans- 

to 



where in (III-5) , the nonholonomic velocities vector, {u j , 

contains the three projections ( o , cd , <0 ) of the angular 

velocity vector onto the body fixed^axil and the three 

projections of the reference point translational velocity 

(u, V, w), onto the body axes. The elements of 7 
k 

(i, j = 1, 2, 3) are direction cosines; the sub-matrix ([y3 is an 
orthonormal rotation transformation relating the attitude of the 
body fixed axis system to the inertial frame. The sub-matrix, 
[^nj, is also a rotation transformation; however, it is not 
orthonormal since it relates vector components based on an 
orthogonal basis to those of a skew basis; namely, the axes 
about which Euler rotations are measured. 


The mass matrix for body k, appears as 




The force/ torque vector, | G | , contains the external forces 

and torques plus any stiffness and damping force that may arise 
through connections with the other bodies making up the system. 
Clearly, | Cj is the implement through which additional system 
forces /torques are "fed back" to the dynamical system. 
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The constraint equations (third of III-2) are written in 
terms of the nonholomic velocities, | u [ , which in turn are 
expressed in terms of the generalized velocities as in (III-5) . 
The coefficients [”] are obtained from expressions of kinematic 
constraint and these same ^ ^ J coefficients are transposed to 
premultiply the vector |X | , providing constraint forces and 
torques . 


B . KINEMATIC COEFFTCIENTS 

This subsection discusses the aforementioned kinematical 
relations involving expressions of relative and absolute 
velocities which lead to the form of the b J coef f icients . The 
discussion will focus on two adjacent interconnected bodies. 



The origins of the body reference systems are labeled m and 
n. The portions of the body where the bodies connect are located 
and labeled p and q. 

In general, for each interconnected pair of bodies, there 
will be five (5) axis systems or coordinate bases. First, there 
will be an axis system fixed to each of the points, m, p, q, n 
(see Figure III-l) . The fifth axis system is a skew or non- 
orthogonal basis comprising direction lines or unit vectors about 
which Euler rotations are measured. The Euler rotations are used 
to describe relative attitudes between the p and q frames. 
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At each connection joint, there will be six (6) components 
of relative velocity (three relative Euler angle rates and three 
relative translational rates that are measured along the skew 
axes) to be expressed. 
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In (III-8), is a (3 X 3) rotation transformation relating 
vector components in the p system to components in the q sys- 
tem. It transforms from p to q. The transformation is 
similar, and note that the^product (R^ R^ = R^) transforms from 
m to q. The matrix, [n]-‘ relates victor components referred 
to orthogonal axes to those referred to skew axes. The matrix 
fs J is a (3 X 3) skew symmetric matrix containing components 
of. m vector positioning point p from m. 
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Finally, it is pointed out that certain rows from (III-8) con- 
stitute rows of [b] and other rows of (III-8) are rows of 
matrix [B] . 


C. INTERCONNECTION CONSTRAINT FORCES 

The interconnection constraint forces/torques can be readily 
obtained by manipulations of the first and last portions of 
(III-2) . Taking the first time deviative of part C of (III-2) 
gives 
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[b] I u ( = I V [ - [b] 1 u I (III-IO) 

using (III-3) in a of (III-2) and together with (III-IO) pro- 
vides an expression for the constraint force/torque vector,^ , 

(III-ll) 

Where we further note certain functional dependences. 

b = bO) 

B = B(0) 

/3 = 0 ( 0 , u) 

S = 5 ( 0 , $ , S , t) 

G = G(u,5 , t) (III-12) 

h = h(0 ,0) 

A = \(0, uj ,6 , t) 
u = u( /3 , u, 0 , 5 , t) 

The variables denoted by 5 and their time deviative, 6 , 

represent control system type variables or other miscellaneous 
variables required to complete the simulation. 

D. ARRANGEMENT OF THE STATE VECTOR 

The variables for the configuration space, y , are arranged 
in vector form and appear as, 
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NB 




P 


N/3 


"N6 
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Where NB represents the number of total ’’bodies" in the system, 
N/3 is the total number of position coordinates necessary to the 
system, and N6 is the total number of auxilliary (control) 
differential equations required. 


Now, given that the | y | vector is known (numerically) 
from^prescribed initial conditions or from numerical integration 
of ^ y [ , the primary task of the solution process is to numer- 
ically establish the | Y [ vector. The | y [vector is numerically 
(step by step) integrated so as to produce an incremented | Y } 
vector, thus a sequence of time point solutions. 


In way of summary, a narrative description of the steps 
(numerical evaluations) necessary to produce | y [ given | Y [ j 
follows . 


The matrices [b J . and b J . are kinematic coefficients 
that depend on position'^and modal displacement variables and are 
evaluated as the first step. 
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Now, if available numerical techniques (also computer 
software and hardware) were absolutely accurate, we would be 
assured that the |u [ . vectors, resulting from numerical integra- 
tion of the { U vectors, would satisfy the constraint Equation, 
c of III-2. Thi^ is not the case; therefore, the second step of 
the solution process is to calculate the dependent elements of 
the I U j- . vectors by using Equation c of III-2. In fact, due to 
anticipating numerical inaccuracies, only the independent element 
of the I U } . vectors are obtained by numerical integration. 

There are otly n-m ^integrators" involved in the solution process 
even though all of the elements of the | U j . vectors are numeric- 
ally evaluated (by use of Equation III-l) ; t?e have good numerical 
resolution in the independent | U } . elements due to using the 
Lagrange multipliers | A ^ 

A kinematically consistent system results from satisfying 
Equation c of III-2. The | U | . vectors may now be used with the 
selection and kinematic transfo^-mations as indicated and Equa- 
tion b of III-2 to produce (numerically) position coordinate 
rates] ft | completing the third step of the process. 


Sufficient calculation has been completed to this point to 
then evaluate the control variable rates as per Equation III-12, 
producing] ^ During the process of calculating the | 5 ] vec- 

tor, all of the required control actuator torques (or forces) 
are calculated, because sufficient numerical information is 


available. All of the constituents of the torques/ force vectors, 
I G j- . , are now available and therefore | ^ » f 

[t] ] are numerically evaluated, (refer to ti\e functional 
expressions of Equation III-12) , x^/hich completes the fourth step 
of the process. With reference to Equation III-ll, we note that 
there is now sufficient numerical information to evaluate | ^ > 

xvhich is then used in Equation III-2 to calculate the { U I , 
completing the fifth and final step of the process. ^ 


It is noted in the above discussions that the solution 
process may be carried out through completion, providing the 
state vector is numerically knox^ml. At any step of a simulation, 
^Le ^y j vector is known, of course, as the result of numerical 
integration. The initial state vector is another matter. It is 
difficult, if not impossible, for a user to prescribe { U }. 
vectors that are kinematically consistent with the conditioAs 
of Equation c of III-2; also, the nonholonomic velocities of 
(U "} . , when considered as a complete set, are of a somex^?hat 
abstract nature. The user ^is in a much better posture to pre- 
scribe initial values of{ /? ) (the initial velocities that are 
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physically meaningful to him). Thus, to initiate the simulation 
(that is, to create an initial state vector from information, 
the user is in a position to prescribe) some preliminary steps 
must be taken as follows. 

The user must prescribe initial values of the {/?},{ P } and 
{ 5 } vectors. Now, in that (the prescribed position rates), 

are explicitly dependent on time and are always available, the 
kinematic equations, b of III-2 and c of III-2, may be used 
together to establish initial values of all i U } 


E. SINGLE TANK/LIQUID MASS COMBINATION 

The specifics that relate to a single tank/liquid mass 
system will be identified in this section. The relevant geometry 
and separate coordinate systems are noted as: 



liquid mass 
triad 


The above sketch depicts the general problem and we will 
make further simplifications which are helpful to our problem. 


For our particular discussions, we will let the body and 
tank axis systems be coincident with each other. Furthermore, it 
will be convenient to introduce a polar spherical coordinate 
system to position the tank mass with respect to the tank axis 
system. Consider the following typical tank axis system. 

Z 


liquid mass 


Y 



X 
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The spherical coordinate transformation is: 

X = r sin cf) cos 0 

y = r sin <|> sin 0 (III-14) 

z = r cos 

The equation for the assumed ellipsoidal surface is, 

("t) +(“t) ' 1 (III-15) 


and for the axi" symmetric case, let b = a; thus the surface con- 
straint equation becomes (in spherical coordinates) , 

a 

— (III-16) 

. 2 . , / a \ ^ 2^ 

Sin 4> + j cos 4> 

and upon taking the first time deviative 




or more concisely, f = <^ ( 4> ) 


(III-17) 


(III-18) 


where it is noted that the condition of axi-symmetry is responsi- 
ble for the dependency of f on and alone. The expression 
given as (III-18) in essence becomes the « expression (III-2c) . 
The following sketch portrays the actual mechanization of our 
problem. 
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ZI 



Note: ^ (eqt. 111-16)= tt - /0V 


ellipsoidal constraint 
surface 


(liquid mass) Body III 


(dummy mass) Body II 



\ ^Body I (tank/ structural mass) 

Body I (frame P 2 ) 


Body II (frame q^) 


-Body III (frame q^) 


Body II (frame p^) 


The frame for the fluid mass (q 3 ) is only permitted to move 
along the radius connecting body I to body III. Body III is a 
ficticious or dummy body that aids in the mechanization process. 
The radius between body I (the tank axis centroid) and body III 
remains constant in length; hence as the "pendulum" rotates the 
fluid mass moves inward and outward in such a manner as to main- 
tain the necessary constraint (ellipsoidal constraint surface). 

It is further noted that a proper choice of Euler permutation for 
positioning frame q^ with respect to frame p^ will lead to an 
explicit expression for both polar spherical angle (and rate) as 
a function of ft and ft within the general dynamical equations of 
motion (III-2) . 

Next, with regard to the expression required for evaluating 
A , given as iii-ii, it is noted that we need the derivative of 
a (givenas finIII-18). To accomplish this, we will use the 
analytical expression for an analog differentiation circuit. 
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K 


1/s 


P 


FIGURE III-2 DIFFERENTIATION CIRCUIT 


With reference to Figure III-2, we observe that, 

P = K( a — P ) 

s 


or 


hence, a = sa== 

Where, for K large, (with respect to s = j<o) , the desired 
expression for a* is identified as the output variable, p , 
from the differentiation circuit. 

Finally, it is pointed out that the tank fluid force is 
obtained from the A vector as that constraint force necessary 
to enforce zero relative translation between the frames p and 
q^. Actually the output from the digital program gives this 
constraint force in the p^ frame. 

The mechanical analog discussed in this chapter has been 
implemented into a digital computer program (LAMPS3) . A users 
guide for the program, with sample input and output, is presented 
in Appendix A, 
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IV, TEST/ANALYTICAL CORRELATION 


This chapter presents a discussion of the Phase I and II 
test results, concentrating on Phase II results. The Interim 
Report (Reference 1) presented Phase I test results. In addition, 
this chapter presents a comparison of test and analytical results 
and comparisons of baffeled and smooth tank data. Appendix B 
presents all Phase II measured force time histories including a 
test log delineating applied accelerations for each test case. 


A. OBSERVATIONS ON LIQUID MOTION 

In addition to providing the force and center of mass motion 
needed for the analytical correlation, the experimental program 
provided some insight into the fluid dynamics of propellant re- 
orientation. A brief survey of the current state-of-the-art in 
propellant reorientation will show how these tests extend the 
technology of this form of liquid motion. 

Propellant reorientation has been studied for many years, 
primarily through the use of drop towers to simulate the low-g 
environment. Various parameters that influence the reorienta- 
tion have been evaluated. Masica and Petrash (Reference 4) and 
Bo\^man (Reference 5) investigated reorientation in a flat-ended 
cylindrical tank while Salzman and Masica (References 6 and 7) 
considered a cylindrical tank with hemispherical domes, one of 
which was inverted. Reorientation in a spherical tank was 
investigated by Labus and Masica (Reference 8). All of this 
experimental work was performed using purely axial accelerations. 
Masica (Reference 9) also investigated the liquid motion pro- 
duced by a purely lateral acceleration. By inclining the tank 
with respect to the acceleration vector, Bowman (Reference 10) 
considered the effect of an off-axis acceleration. However, the 
applied acceleration was still perpendicular to the initial gas/ 
liquid interface. The initial condition for all of Bowman’s 
tests was a flat gas/liquid interface, while in the other tests 
mentioned here a highly curved, near zero-g interface was 
established before applying the acceleration. 

In all of these investigations the Bond number was used to 
categorize the liquid reorientation. It is the ratio of the 
gravity force to the surface tension force. 
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where: P = liquid density 

a = acceleration 
r = tank radius 
^ - liquid surface tension 

When the interface is initially flat and the Bond number for 
the reorientation is less than 10, the motion of the liquid is 
along the tank wall. At larger values of Bond number, a central 
instability forms, whose size is a function of the Bond number. 
With this condition, a large portion of the liquid reorients 
through the center of the t^nk. However, with the initial low-g, 
highly-curved interface, the liquid motion was always totally 
along the tank wall for the range of Bond numbers evaluated 
(3 to 450). 

When the tank was inclined at angles as small as one degree, 
the geometry of the tank relative to the interface influenced 
the liquid motion. The instability still formed in the center 
of the interface, but was laterally displaced so most of the 
flow was eventually along one side of the tank. Apparently the 
difference in interface curvature on opposite sides of the 
instability caused the force that produced its lateral displace- 
ment. 

Another noteworthy experimental program was performed by 
arresting a free- falling tank (Reference 11). This allowed a 
larger tank and larger accelerations, yielding Bond numbers as 
large as 5800. When the acceleration was axial, the liquid motion 
was in the form of rain, traveling through the center of the 
tank. By laterally oscillating the liquid immediately before 
the test, off-axis effects were considered to some extent. At 
times, the liquid motion was along one side of the tank, depending 
upon when the tank was arrested in the slosh cycle. 

The effect of ring baffles on the axisymmetric reorientation 
of the liquid in a cylindrical tank has also been evaluated 
(Reference 4) . 
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1* Test with Baseline Tank - All of the Phase I tests, in 
this study, were conducted with the large tank which has smooth 
internal walls. In all cases, a lateral acceleration, acting 
parallel to the initial liquid interface, was applied. The 
test provided an evaluation of the influence of a true off-axis 
acceleration on the manner of reorientation. 

In every case the liquid reoriented in an un symmetrical 
manner, along one side of the tank. This was true for a wide 
range of lateral accelerations. In some of the tests, the lateral 
component was estimated to be 0.005g and lasted for a period of 
only 0.3 second at the beginning of the test. Compared to the 
axial component, the lateral component ranged from 5 % to 50%. 
Regardless, the manner of liquid motion remained essentially 
the same for all Phase I tests. 

In those tests with the larger axial acceleration, the initial 
formation of an instability in the center of the interface 
was observed. It appeared as a hump in the surface that gradu- 
ally joined the wall flow and disappeared. The size achieved 
by the instability was dependent on how long it had to grow 
before joining the wall flow. 

As the liquid began to move, the liquid interface remained 
relatively flat so the motion appeared as a rotation of the 
interface about its center. This was most pronounced at the 25% 
and 50 /o liquid volumes. Once the leading edge of the flow 
reached the tank dome, the liquid interface began to acquire 
some curvature. In general, the liquid overshot its final 
equilibrium position, continuing on around the tank and recircu- 
lating a small percentage of the liquid. The liquid center of 
mass traveled beyond the equilibrium position and was coming to 
a stop as the test ended. If the test could have been continued, 
damped oscillation of the liquid about its equilibrium position 
would have been observed. 

Very little splashing of the liquid was observed during the 
reor ien tat ion. The leading edge of the liquid strongly adhered 
to the tank wall and any turbulance was confined to the surface. 

A typical test is shown in Figure IV-1. A small instability 
hump can be observed forming in the third photo, but by photo 4 
it has joined the wall flow. In photo 5, the liquid has spread 
out and some of it is recirculating back to the initial position. 
At the end of the test (photo 6), the liquid is beyond the equili- 
brium position and some liquid is still being recirculated. 



IV- 4 


In all of the photos in this section, the axial acceleration 
is acting downward and the lateral acceleration acts to the left. 

When the liquid volume was 75%, the reorientation was similar 
to that described above, except that the ullage assumed the form 
of a bubble and moved to the opposite end of the tank. The bubble 
followed the tank wall as it moved. Its surface was highly 
irregular due to the flow of liquid around it. At the end of 
the test , the bubble had become somewhat flattened and the center 
of gravity of the liquid had overshot the equilibrium position. 

During the Phase II test program, two tests were performed 
in which only an axial acceleration was applied to see if the 
reorientation would be symmetrical. In one of the tests, the 
acceleration was in alignment with the tank axis. A sequence of 
photos from this test are shown in Figure IV-2. As expected, an 
instability formed in the center of the interface that had the 
shape of a hollow cylinder. By photo 3, the instability had 
become unsymmetrical and by photo 5, the bubble within the 
instability had gro\^ so that most of the liquid was on the walls 
of the tank. At the end of the test (photo 6), the liquid was 
collected but was turbulent and numerous gas bubbles had been 
entrained. Some of the liquid was returning to the tank bottom. 

In the second test, the tank was inclined at 45^. The motion 
was similar to the above test but there was more distortion of 
the instability and more liquid recirculation due to the geometric 
configuration. 

Based on these tests, it can be concluded that symmetric 
reorientation with its attendent instabilities and geysers as 
have been observed in past test programs, is a special case that 
is not very likely to occur in practice. Any lack of symmetry, 
caused by geometry or acceleration, will cause unsymmetrical 
reorientation of the liquid. 

2. Tests with Baseline Tank Containing Baffles - Adding ring 
baffles to the large tank produces liquid motion much different 
than that observed in the tank with smooth internal walls. Three 
examples of the influence of the baffles are shown in Figures 
IV- 3, IV-4 and IV- 5. 

As the leading edge of the liquid comes into contact with a 
baffle, the liquid is strongly deflected. Photos 2 and 3 of 
Figure IV-3 are the best example. The liquid was deflected by 
the lower baffle and it passed over the other baffles without 
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touching them. Photo 5 of the same test shows another deflection 
of the liquid after it had crossed the top of the tank. Some 
liquid was deflected back toward the tank bottom, but most of it 
was confined at the top of the tank by the upper baffle. 

For the test shown in Figure IV-4 the leading edge of the 
liquid starts at the upper baffle but the liquid was still de- 
flected away from the tank wall (Photo 2). At the end of the 
test there was more recirculation due to the relative orienta- 
tion of the tank and the acceleration. Some liquid passed over 
the baffles and returned to the bottom of the tank (Photos 4 
and 5) . 

In Figure IV-5 there was little deflection of the liquid 
as it first began to move. There was some deflection of the 
liquid (Photo 5) as it hit the baffles on the other side of 
the tank. However, this deflected liquid quickly recirculated 
so that most of the liquid was collected by the end of the 
test (Photo 6) . The baffles were the most effective in this 
case and a minimum of turbulence was induced. 

The motion presented in these photos is typical of the 
motion observed in all the baffled tank tests. As the leading 
edge of the liquid impacts a baffle, a significant deflection 
occurs. The baffles act to reduce the recirculation of the 
liquid that was observed in the tests with the smooth walled 
tank. It appears that the collection of the liquid at the 
equilibrium position is speeded, but added turbulence is in- 
duced . 

3. Tests With Half-Sized Tank 

The effective test time was increased by using the small 
tank which is one-half the size of the large tank. Two of the 
small tank tests are shown in Figures IV-6 and IV-7. 

With the large tank, the liquid was observed to reorient, 
overshoot its equilibrium position and begin to come to rest. 

With the small tank the liquid reoriented, overshot the equili- 
brium position and came to a stop, and then moved back towards 
equilibrium again. It was nearly at rest at the equilibrium 
locationat the end of the test. It appears that if additional 
test time were available there would be a small amount of liquid 
motion as it stabilized in its final position. 
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The portions of the reorientation that are common between 
the large and small tank appear the same. For example, compare 
the photos of Figure IV- 1 with the first five photos of Figure 
IV-7. In both Figures IV-6 and IV-7, Photo 5 shows the liquid 
as it came to rest after overshooting the equilibrium position. 
Photo 6 shows the liquid as it nears the equilibrium position 
at the end of the test. 

In all the tests the liquid became well collected as it 
was returning to the equilibrium position following the over- 
shoot. Any liquid that broke away from the leading edge and 
recirculated, continued on around the tank and rejoined to form 
a single liquid mass . 



REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 






IV-8 



I. 0. 0 sec 



3. 0. 6 sec 




2. 0.4 sec 



4. 0. 8 sec 



6. 1. 6 sec 


Rgure IV-2. Test I , Phase II ; 25 % fill ; ex= 0° ; Aa= 0.084 ; no lateral acceleration 
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I. 0.0 sec 2. 0.4 sec 



5. 1.0 sec 6. 1.7 sec 




Figure IV-3. Test 6 , Phase II ; 25 % fill ; ex= 0° ; Aa= 0.074 
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5. 1. 7 sec 


Figure lV-4. Test 7 , Phase II ; 25 % fill ; qx= 45° ; Aa= 0.074 
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Figure IV-5. Test 8 , Phase II ; 25 % fill ; ex= 90 ; Aa=0,075 









Test 15 - Phase II , 25 % fill 
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B. DISCUSSION OF TEST/ANALYTICAL CORRELATION 


This section presents a qualitative discussion of several 
aspects of the analytical model based on analysis of the test 
force data and photographic records. In particular, the con- 
cepts of effective liquid mass and fluid damping will be dis- 
cussed in detail. Subsequent sections will present comparisons 
between analysis, test, Phase I and Phase II results. 

1. Liquid Effective Mass - Phase I test/analytical corre- 
lation, as presented in Reference 1, always showed analytical 
force predictions to be of a larger magnitude than measured 
forces . The force exerted on the tank by the moving liquid is 
composed of two components: 1) a centripetal acceleration 

force due to liquid velocity, and 2) the D’Alembert force due 
to the applied acceleration field. In addition to this inertial 
reactive type force, there may also be a viscous dissipative 
force due to friction between the moving liquid and the tank 
wall. 


During Phase I investigations, it was thought that the 
overall magnitude of liquid forces could be reduced by increas- 
ing this viscous dissipative force, thereby reducing the liquid 
velocity and .decreasing the cenLripetal acceleration component 
of force. Subsequent analyses indicated that, indeed, the 
liquid forces were reduced, however, the decrease in liquid 
velocity degraded the time correlation of peak forces when 
comparing analytical and measured data. Increasing the viscous 
dissipative force enough to reduce the analytical force peaks 
to match measured data resulted in unacceptable time delays in 
the peak analytical force. Initial Phase II studies attempted 
to define this viscous dissipative force in a manner which 
would improve correlation in force magnitudes while maintain- 
ing time correlation. These efforts proved unsuccessful. 

Damping investigations eveiiLually lead to the development 
of the concept of liquid effective mass. Indepth studies of 
photographic records of the drop tests revealed that during 
liquid reorienta cion , the liquid begins to assume a somewhat 
curved interface as it accelerates, resulting from expansion 
of the liquids surface area. This characteristic can be ob- 
served in Figures IV- 1 , IV-6 and IV-7. Analytically breaking 
the liquid into finite sections, as shown in Figure IV-8, re- 
veals Lhat eacii segment exerts a centripetal acceleration force 
on the t£mk , perpendicular to its velocity vector. From Figure 
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IV-S, it is obvious that the net centripetal acceleration force 
exerted on the tank is less than that obtained by lumping the 
total liquid mass at the liquid center of mass (cm). 



F. 

1 


< F 


direction 


cm 


where : 


= unit normal vector to 
the constraint surface 
at the liquid cm 

F^ = segment centripetal ac- 
celeration forces 

F = centripetal acceleration 
cm ^ 1 . 

force resulting from 

lumping total mass at 

liquid cm 


Figure IV-8. Liquid Effective Mass Justification 


Based on test results, the effective mass appears to be a func- 
tion of tank percent fill volume. It may also be a function of 
tank L/d, however, all tests were run with similar aspect ratio 
tanks so this functionality cannot be assessed here. 

Figure IV-9 delineates the liquid effective mass factor 
which when applied to the centripetal acceleration component 
of liquid force, best reproduces the test results. The effec- 
tive mass factor (expressed as a percentage of the total liquid 
mass) is presented as a function of tank fill volume. 
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% TANK FILL VOLUME 


Figure IV-9. Liquid Effective Mass Factor 


The moving liquid force can be expressed as: 


F = MASFAC * F 


CENT 


+ F 


D'ALEMBERT 


(IV- 1) 


where: MASFAC = effective mass factor (Figure IV-9) 

= centripetal acceleration force component 

CENi 

^D'ALEMBERT = D'Alembert force component 


Note that the effective mass factor is only applied to the cen- 
tripetal acceleration component of liquid force; the velocity 
dependent part. When the liquid totally reorients and comes to 
rest, velocity equals zero, the force exerted by the liquid is 
equal to: 


\lQUID ^APPLIED ’ 


(IV-2) 


the D'Alembert force, where H^XQUID total liquid mass. 

Application of the effective mass factor in the analytic 
model resulted in good correlation between analytical and 
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measured force magnitudes, while maintaining good time corre- 
lation of peak forces. The curve of effective mass factor in 
Figure IV-9 is based on tests ranging from 10% to 75% tank fill 
volumes. As the fill volume increases, the centripetal accel- 
eration force component decreases drastically due to the nature 
of the reorientation. At 75% fill volumes, the liquid was ob- 
served to flow around a bubble formed by the ullage gas which 
tended to traverse the tank wall. This behavior results in 
smaller liquid velocities and, therefore, decreased centripetal 
forces. On the other hand, in the 10% fill case, the centri- 
petal force tended to approach the value based on the total 
liquid mass. Intuitively, based on test results, it is doubt- 
ful that an effective mass factor of 100% is ever reached. 

2. Liquid Dam p in o; - Two criteria exist for evaluating the 
worth of the analytical model; one is the force time history 
comparison and the second is how well the analytic liquid cm 
position compares with the liquid motion photographically re- 
corded during the tests. As discussed above, the effective 
mass factor provides good force magnitude correlation, and 
peak force timing appears to correlate well without any viscous 
dissipative forces applied. However, a study of the drop test 
films indicates that, through some mechanism, the energy of the 
moving liquid is dissipated. It appears that the liquid is 
slowed by a continuity phenomena. That is to say, as the 
liquid velocity builds, the liquid crossec tiona 1 area decreases, 
increasing the liquids surface area. This expansion of the 
liquids* surface area appears to convert flow energy into 
strain energy (surface area tends to maximize), thus slowing 
the liquid. This expansion is especially evident as the liquid 
traverses the tank dome. The expansion appears to include some 
out of plane flow since the liquid in the tests was not physi- 
cally constrained to two dimensions. As the liquid slows, a 
contraction process occurs, resulting in minimum surface area 
as the liquid comes to rest. This expansion and contraction 
is best delineated by the small tank reorientations shown in 
Figures IV-6 and IV-7. 

Phase II comparisons of liquid cm position compared to 
film records of various tests, indicated that some dissipative 
force was required to maintain a reasonable liquid cm time his- 
tory correlation. Obviously, a point mass model cannot hope to 
actually represent the expansion and contraction phenomena dis- 
cussed above. However, a compromise has been found which 
appears to maintain a realistic center of mass location, during 
reorientation, without degrading force magnitude or force/ time 
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correlation. The photographic records indicate that the liquid 
expands and slows the most as it traverses the tank dome. To 
represent this phenomena, a dissipative force of the form, 
f = /, has been incorporated in the analytical model. In 

this formulation j8 is the angular rate (RADIANS/ SEC) of the 
liquid cm as it orbits the center of the tank on the ellipti- 
cal constraint surface (see Figure IV- 10). The parameter, M, 
is the **damping" coefficient. It has been found that a value 
of M = .00002, provides satisfactory center of mass location 
correlation for both the baseline tank (Phase I) and the half 
size tank (Phase II). This correlation is presented in Figures 
IV- 11 through IV- 16. In these tests the analytic liquid cm 
position maintains fairly good correlation with the intuitive 
cm location. Since the dissipative force, f, is not intended 
to represent a wall friction effect, it is not included in 
liquid force calculations. It is only used as a force on the 
liquid cm to dissipate flow velocity. It is felt that liquid/ 
wall friction forces are extremely small compared to inertial 
reactive forces between the tank and liquid. 



liquid dissipative 
force , f = /^/ ^ / 

inertial reactive 
force between the 
tank and liquid 

liquid velocity 
vector 


Figure IV- 10 . Liquid Damping Visualization 


It is not known whether the value of ^ used in analysis 
of the test cases is applicable for a large class of tanks. It 
would appear that the value would have to be a function of tank 
size and perhaps other variables. However, based on the limited 
parametric variations in the Phase I and Phase II tests and the 
small size of the tanks, a definitive evaluation of the func- 
tional characteristics of M could not be conducted. Further 
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data utilizing much larger tanks is required. However, it 
should be noted that neglecting the dissipative force, all 
together, results in conservative force magnitudes which may 
be sufficient for some analyses . 

3. Constraint Surface Definition - The mechanical analog 
developed during Phase I investigations portrayed the liquid 
as a point mass moving on a constraint surface which was deter- 
mined by rotating the tank (analytically) in a one-g field; the 
constraint surface was defined as the locus of center of mass 
locations prescribed during the rotation, assuming the free 
surface was planar. This constraint surface was approximated 
by piecewise continuous elliptical segments. When the liquid 
cm deviated some tolerance from the constraint surface, the 
elliptical coefficients were updated to return the liquid cm 
to the desired surface. During Phase II it was decided to 
alter the constraint surface definition slightly. It was found 
that updating the elliptical coefficients resulted in slight 
analytical force discontinuities. To correct this problem, 
the constraint surface is now defined as the ellipsoid which 
best fits the constraint surface discussed above. Hence, for 
any given tank geometry and fill volume, the constraint surface 
is approximated by a single set of elliptical surface coeffi- 
cients; no updating is performed. Comparisons have shown that 
this does not substantially alter the analytical liquid forces, 
but does clean up the analytical results. 

Test/Analytical Comparisons - Bare Tank - Utilizing 
the concepts of effective liquid mass and liquid damping dis- 
cussed above, the analytical model (computer program LAMPS, 
Reference 1) was used to simulate several Phase I (large tank) 
and Phase II (small tank) tests to assess the degree of corre- 
lation. Figures IV- 11 through IV- 13 present the correlation 
for three Phase I (1974) tests: 13, 16 and 17. Tests 13 and 

16 (Figures IV- 11 and IV- 12) were conducted using the baseline 
tank (2.5 inch radius) and a 257o fill volume. Test 17 (Figure 
IV- 13) was conducted with the same tank and a 50% fill volume. 
The figures present an overlay of analytical and test force 
time histories in the Y and Z directions (tank axis system) . 

In addition, sketches are presented which depict the liquids* 
cross sectional shape, traced from film recordings, overlayed 
with the analytical liquid cm position. For these tests the 
correlation of both force and liquid cm position is very good. 

In the Phase II test program, a half sized tank (1.25 
inch radius) was used to repeat several of the 1974, Phase I 
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tests. Use of the half sized tank results in a greater equiva- 
lent test time. Figures IV-14 through IV- 16 present analytical 
correlations for three of these small tank tests. Tests 15 and 
16 (Figures IV-14 and IV- 15) were run using a 257o fill volume 
and correspond to tests 13 and 16 (respectively) in the Phase 
I test program using the large tank. Test 21 (Figure IV- 16) 
was run with a 257o fill volume and approximately twice the 
lateral acceleration of Test 16 (Phase II). 

The forces due to the moving liquid in the small tank 
tests are much smaller than those in the Phase I tests. Con- 
sequently, the sensitivity settings of the charge amplifiers 
for the load cells were greatly increased to insure adequate 
force resolution. The load cells are susceptible to thermal 
drift from the drop capsule lighting system and the increased 
sensitivity, required for the small tank appears to have re- 
sulted in some drift, even though the load cells were insulated. 
Figures IV-14 through IV- 16 show fair force correlation for the 
small tank tests. 
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However, the cm location timing appears to be excellent. It is 
felt that the small force magnitude in these tests coupled with 
thermal drift tend to make the measured data suspect. However, 
based on the large tank tests and the center of mass correla- 
tion on the small tank tests, it appears that the analytical model 
provides reasonable results even for the small tank. Study of 
Test 21 results (Figure IV- 16) indicates similar force trends 
between measured and predicted forces. This test of course em- 
ployed larger lateral acceleration's and, hence, larger liquid 
forces . 

Overall applicability of the analytical model concept of 
a point mass moving on an ellipsoidal constraint surface is 
considered to be excellent. The test correlations with the 
bare tanks (Phase I and Phase II) indicate that the analytical 
model provides a reasonable approach for predicting forces due 
to transient liquid motion. 

5. Discussion of Baffled Tank Test Results - During Phase 
II testing twelve drops were performed using the Phase I tank 
(2.5 inch radius) with three ring baffles installed (see Chapter 
II) . It was initially hoped that baffles could be analytically 
represented by some form of drag or friction force within the 
framework of the current point mass model. However, the ob- 
served liquid motion during the Phase II tests indicated that 
this approach was not feasible. During the low g reorienta- 
tions simulated, the baffles, when encountered, essentially 
redirected the liquid flow to the tank interior. Figures IV- 17 
through IV-21 present the results of five of the baffled tank 
tests. From the sketches of liquid motion, it is noted that 
some liquid is trapped under the baffles and laterally tra- 
verses the tank. Eventually, at second contact with the baf- 
fles, the liquid becomes a spray. The simulation of these 
complex flow patterns is not feasible with a point mass model. 




















Figure IV-21. Baffled Tank Comparison 
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Although the analytical model cannot simulate the flow for 
the baffled test cases, much was learned about liquid reorien- 
tation from these tests. A study of the force data from these 
tests has also yielded some interesting observations. Many of 
the baffled test cases repeated the test conditions used in the 
unbaffled tests of Phase I. In Figures IV- 17 through IV-21, in 
addition to liquid motion sketches, force time history compari- 
sons are shown between the particular baffled test case and its 
corresponding Phase I, unbaffled test case. It is interesting 
to note the extreme similarity in measured data. In general, 
the baffles tend to somewhat reduce the force levels and delay 
the peak forces. However, this is only a general qualitative 
evaluation. It appears, however, that using the analytical 
model and neglecting the baffles, would result in fairly rea- 
sonable force predictions. The results would at least be con- 
servative in most cases. Larger fill volume cases (i.e., 50%) 
exhibit flow which better approximates the flow in unbaffled 
cases. In these cases, the baffles might be represented as an 
increase in liquid damping force, at least during the initial 
part of the reorientation. 

6 . S ymmetric Reorientation - The Phase II test program 
included two tests using the Phase I tank (2.5 inch radius, no 
baffles) in which only axial acceleration was applied. Section 
A of this chapter discussed the observed reorientation in these 
cases. It has been found, in this study and others, that sym- 
metric reorientation, where the liquid does not follow the tank 
wall, but traverses the tank interior, is a special case which 
will not generally occur in practice. This type of reorienta- 
tion requires that the liquid is completely at rest and the 
applied acceleration must be perpendicular to the liquid sur- 
face and along the tank axis of symmetry. Test 1 of the Phase 
II test program satisfied these requirements. The reorientation 
is shown by the photographs in Figure IV-2. Due to the limited 
occurance of this type of reorientation and the limited test 
data available for correlation, the analytical model was not 
structured to simulate these cases. However, the measured 
force data does provide some useful information. 


Previous studies have been conducted on impact forces on 
model propellant tanks produced by symmetric reorientation 
(Stephens, Reference 11) . These tests were conducted at higher 
acceleration levels than those used in this program. Stephens 
found that impact forces on the order of 2 to 3 times the hy- 
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acceleration on the order of 0.2 to 0.5 g. For Test 1 of this 
study, however, it was found that the measured force was approxi- 
mately equal to the hydrostatic force. During reorientation 
the axial force on the tank slowly increased as the liquid 
traversed the tank or collected at the tank dome (Figure IV- 2) ; 
no impact type force was observed. Figure IV-22 delineates 
these results for Test 1. 

Based on the contradiction between the results of Refer- 
ence 11, and Test 1 of Phase II, it appears that more tests 
would be helpful in understanding symmetric reorientation in 
low g environments. However, since this type of reorientation 
is considered a special case, resolution of the problem is 
mostly of academic interest. 


FZ (Pounds) 



Time (Seconds) 


Figure IV-22. Symmetric Reorientation Forces 
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V, CONCLUSIONS AND RECOMMENDATIONS 


Transient liquid motion, in low g environments can be 
accurately simulated by an analytical model which represents 
the liquid as a point mass moving on an e llipsoida 1 cons traint 
surface. The concepts of liquid effective mass and liquid 
damping developed in the Phase II study have greatly improved 
the accujracy of the model*s force predictions and center of 
mass location cor re la t ion. 

The mechanical analog has been cast in a framework which 
allows the incorporation of tank/liquid systems into general 
spacecraft simulation problems. The three dimensional model 
has been structured to facilitate incorporation into general 
spacecraft equations for use in a variety of analyses: i.e., 

control interaction during orbital maneuvers, reentry analy- 
ses, loads analyses, etc. 

1. The mechanical analog should be integrated into on- 

going and future analytical efforts to improve the 
simulation of liquid/structure interactions: i.e.. 

Shuttle ET reentry studies. Space Tug control sys- 
tem design and maneuver analyses, etc. 

The test system developed during this program is capable 
of providing insight to the character of liquid reorientation 
and the forces exerted on Spacecraft by moving liquid. It 
provides the capability to apply arbitrary low g accelerations 
while recording transient liquid forces and providing photo- 
graphic documentation of the test. 

2. Further testing should be conducted to build the data 
bank necessary for further analytical model verifica- 
tion. Testing should include various tank geometries; 
ogive, conical, etc. Additional scale propellant 
management devices should be incorporated into the 
tests ; various baffle arrangements , surface tension 
management devices, etc . 

Reorientation in baffled tanks has been found to be a 
complicated phenomena. Simulation of the observed flow patterns 
with a point mass model is not feasible. However, comparison 
of similar tests with and without baffles revealed a striking 
similarity in measured force trends. Baffles, in general, tend 
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to reduce force levels and delay (time wise) force peaks. 
Based on these observations, it appears that the analytical 
model may be used, neglecting the baffles, for conservative 
force predictions in some applications. 
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APPENDIX A - THREE DIMENSIONAL MECHANICAL ANALOG COMPUTER PRO- 
GRAM USERS GUIDE (LAMPS3) 


« This appendix discusses, in detail, the implementation of 

the tank/liquid system into the framework of the DYNAMO com- 
puter program (Reference 3). Chapter III discussed the general 
analytical techniques employed in DYNAMO: in this appendix we 

shall concentrate on the specifics of the mechanization, in 
particular the required input data. 

Method of Mechanization 


The implementation of the mechanical analog into DYNAMO 
required some significant alterations to the basic computer 
code. In particular, the specification of a state dependent 
constraint surface (ellipsoidal) required the addition of a 
capability to differentiate a (see Chapter III, Section E) . 
The modified DYNAMO program has been named LAMPS3 (Large AMP li- 
tude ^losh _3 dimensional). This version only contains the non- 
linear time response options in DYNAMO and is restricted to an 
interconnected system of rigid bodies . 

The specific implementation of the tank/liquid mechanical 
analog is shown in Figure A-1. 



z 
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\ 



* 


1, Drop Capsule 
Mass 


Body 3, Liquid Mass 


Y 

p^ Body 2, Dummy Mass 


FIGUBUE A- 1 . TANK/LIQUID MECHANICAL ANALOG 


In the simulation, Body 1 represents the drop capsule mass 
(ie, rigid spacecraft). Body 2 is a dummy mass required for 
mechanization of the mechanical analog in the framework of 
DYNAMO. Body 2 remains a fixed distance from Body 1 during 
the simulation. Body 3 represents the liquid mass and is con- 
strained to move along the vector from Body 1 to Body 2. Body 
3 slides along this vector in order to remain on the ellip- 
soidal constraint surface. The bodies are connected at "hinge 
points". Bodies 1 and 2 are connected at the hinge represented 
by the origin of coordinate systems P 2 and q 2 (Figure A-1). 

At this hinge, translation between the frames is prohibited 
but there is free rotation. Bodies 2 and 3 are joined by the 
hinge represented by coordinate systems P 3 and qj. Frame q 3 
is only allowed to translate in the Zp^ direction, all other 


EEPRODUdBlLITY OF THE 
ORIGE^AL PAGE IS POOR 


relative motion between the frames is prohibited. Hence, it 
can be seen that the simulation represents a variable length 
pendulum, free to rotate at the hinge P2/q2- The length of 
the pendulum varies such that Body 3 (the fluid mass) always 
remains on the prescribed ellipsoidal constraint surface. Ex- 
ternal forces (Fx, Fy, Fz) are applied to Body 1. These forces 
represent the spring motor forces applied to the drop capsule 
during the test program (see Chapter II). In a more general 
sense, these forces could represent reaction control system 
forces acting on an orbiting spacecraft. 

The position of the liquid mass at any time can be ex- 
pressed in terms of spherical coordinates relative to the tank 
coordinate system (frame P 2 , Figure A-1). Figure A-2 delin- 
eates the spherical coordinate set used in the simulation. 



FIGURE A-2. SPHERICAL COORDINATE SET 


The constraint surface is assumed to be an axisymmetric 
ellipsoid. The equation of this surface can be written: 


2 2 2 
ax + ay +cz =1 


(A-1) 


where a and c are the ellipsoid surface coefficients. 
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Transforming equation (A-1) to spherical coordinates and solving 
for R, yields: 



^ 1/ a 

2 

sin 0 + 

— cos^0 
a 



(A- 2) 


Differentiating equation (A-2) results in the required rate of 
change of length of the pendulum in our simulation. 


- 4 I” sin0 COS0 1 - 

1 a 
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. 2^ c 2^ 

sin 0 H — cos 0 
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3/2 


(A- 3) 


Note that R is only a function of 0 and 0 due to the Z axis 
symmetry of the constraint surface. R is the a discussed in 
Chapter III (equation III-2c) - As pointed out in Chapter III, 
it is necessary to compute V , or R in our case, in order to 
correctly calculate the A*s; forces exerted by the liquid on 
the tank. The analog differentiation circuit discussed in 
Section E of Chapter III was employed for this purpose. That 
is, R was expressed as per equation III-19, 


R = K R (1 + |) 


(A-4) 


For values of K large compared to S = j o , where a> is the 
system response frequency, R simply becomes. 


R 


K R 


(A- 5) 


This algorithm has been found to work exceedingl^^^^well for our 
application. LAMPS3 outputs plots of R(t) and jR(t) for com- 
parison with the explicitly defined R(t), equation (A-3) (see 
the sample output at the end of this appendix) . A value of 
K=200 was found to be adequate for this simulation. Note that 
the plot of R is not absolutely smooth. This slight roughness 
is reflected in the liquid forces. However, the effect of 


this roughness on the overall simulation is negligible and 
greatly outweighed by the increased capability achieved by 
implementing the mechanical analog in the more general frame- 
work of DYNAMO. The value of K=200 can be overridden by the 
input data if required to achieve a better definition of R. 

« • • 

In the context of our simulation, by specifying R and R, 
we are really specifying the motion of the frame relative 
to the P 3 frame in the Zp 3 direction (Figure A-1, Z 3 / 2 ) . Note 
that the rate of change of R and Z 3/2 are just compliments of 
each other. 

The radial motion of Body 3 (R) requires the definition of 
0 and 0, the spherical coordinate angle. The values of 0 and 
0 are not readily available in the mechanization. However, 
the Euler angle rotations and rotational rates of frame q 2 
relative to frame P 2 are available. These are the (i 's and 
P's expressed in equation III-2 (Chapter III). From these 
P 's and p's, 0 and 0 (also 9) can be determined as follows. 

Figure A-3 depicts unit vector triads in the P 2 and q 2 
frames . 



FIGURE A-3. RELATIVE ROTATIONS BETWEEN q 2 AND P 2 TRIADS 


Program DYNAMO allows the user to select the Euler angle permu- 
tation order desired. We have chosen the order: 


The rotation 0 has also been constrained to zero since spin- 
ning about the radius, R, has no meaning. Based on this selec- 
ted order, we can express the unit vectors of the p triad in 
terms of the q triad via the direction cosines. 
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Based on the definition of the angles 0* and 0^^ in Figure A-3, 
we can now write the following equations. 


^ = cos 0* = cos P - cos P ^ a 

q P 13 


K • L = sin cos 0^f = sin P^ cos P^ b 
q p 13 


"K • t ~ sin 0* sin 0* = -sin p c 

q P 3 


Hence, from equation (A-6)a: 


0* = cos ^ (cos P^ cos P , 


(A- 6) 


(A- 7) 


ESL-KODUC’i-n.."^^ OF THE 
OElCmAL PAGE i3 POOR 


and ratioing parts b and c of equation (A-6), 


0* = Tan 


sin 13^ cos p ^ 


(A-8) 


Differentiating equation (A-7) results in. 


0 * = 


sin P^ cos P^ + P^ cos P^ sin P^ 
sin 0* 


(A- 9) 


All of the P ‘s and P *s required to calculate 0*, 0* and 0* 
are known and available within the LAMPS3 program. 

Note that the magnitude of the spherical coordinate rate 
0 (Figure A-2) is the same as 0* (equation (A-9)). The sign 
of 0 is based on whether Body 3 is going toward or away from 
the Zp 2 axis. By observation, we can also define the spheri- 
cal coordinates 0 and 9 as; 

0 = 180° - I 0*1 

and 0 = 0* - n 


0°< 0* < 180° 
0 °< 0 * < 360 ° 


LAMPS3 also calculates these variables (0, 0, 0) and plots 
their time history to facilitate locating the liquid cm at 
any time. 

Liquid damping, as discussed in Chapter IV (Section B2) , 
is easily implemented within LAMPS3. Within the code, an 
allowance is made to specify rotational dampers at the "hinges". 
Liquid^ damping^ is achieved by specifying damping proportional 
to /?i(0y) and /?2_(0x) at hinge p2/q2> Damping (m) is assumed 
proportional to 0, and based on equation A-9, the equivalent 
damping on /j , and /? 2 calculated as follows. 
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M sin 13 ^ cos /? 2 
sin 0 



fjL cos f3 ^ sin ^ 2 
sin 0 


(A- 11) 


The concept of effective mass (discussed in Chapter IV-Bl) 
is implemented as follows. As discussed in Chapter IV, the 
effective mass factor (equation IV- 1) is only to be applied to 
the centripital acceleration component of force. Separation 
of the centripital acceleration component of force from the 
other components by altering the basic kinematic equations of 
the DYNAMO program is not feasible. A simpler method of in- 
corporating the effective mass factor is to apply an external 
force to the liquid mass which cancels the desired proportion 
of the centripital force. Figure A-4 delineates the applica- 
tion of this force. 


Constraint 

Surface 



Q = angular rate, rad/sec 

R = pendulum length 

^EM ” sppli^d force required 
to reduce centripetal 
acceleration force 


FIGURE A-4. LAMPS3 IMPLEMENTATION OF EFFECTIVE MASS FACTOR 


The force Fem in Figure A-4 is defined as: 


^EM ^ 


where: 

MASFAC 


R 

n 


(l-MASFAC) R Q ^ 

= liquid mass 
= effective mass factor 
= pendulum instantaneous length 
= angular rate of the pendulum 


(A- 12) 


This applied force is opposite in direction to the centripital 
acceleration force exerted by the liquid on the tank, and 
hence, cancels the undesired portion of that force. The ap- 
plication of this force is automatically built into LAMPS3. 

Its magnitude is determined by the value of MASFAC supplied 
by the user in the input data. 


Input Data for LAMPS3 

The input data for LAMPS3 follows the basic format for 
DYNAMO (Reference 3, Volume II). Once the input data has been 
established, only a few of the values ever change; the rest of 
the input data always remains the same. Part of the input data 
requires knowledge of the initial position of the liquid cm. 

In order to preserve the input format used in DYNAMO, a separate 
program (TANK) has been coded which generates the initial liquid 
cm position, constraint surface coefficients, and liquid mass 
for input to LAMPS3. Table A-1 delineates the required input 
data for this program. A sample output from program TANK is 
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TABLE A-1. PROGRAM TANK INPUT DATA FORMAT 


I NPUT DATA 
999 


READ(A6,4X, 3A6)RUNNO,UNAME 
IF (RUNNO , EQ . 4H STOP) STOP 


10 


READ(12A6) 

READ(12A6) 

CALL COMENT 
RE AD (5E 10. 3) 

READ (3E 10. 3, 21 5) 

IF(ITRIAD.EQ.l)GO 

READ(5E10.3) 

GO TO 999 


TITLE 1 
TITLE2 

COMMENT CARDS, LAST CARD=10 ZEROS COLS 
XL , TR , TD , PC VOL , FDEN 
G(l) ,G(2) ,G(3) ,ITRIAD,ITYPE 
TO 10 

ANG(l) ,ANG(2) ,ANG(3) 
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DEFINITION OF INPUT VARIABLES 


RUNNO 

UNAME 

TITLE 1 

TITLE2 

XL 

TR 

TD 

PCVOL 

FDEN 

G(l-3) 


ITRIAD 


ITYPE 


Run number printed in page headings 
User name 

Title card printed in page heading 

Title card printed in page heading 

Length of propellant tank cylindrical section 

Tank radius 

Height of tank domes from end of cylindrical section 
Percentage tank fill. LE. 100. 

Liquid mass density 

Initial acceleration field acting on tank to cause 
initial liquid positioning . 

G(l) = X axis acceleration 
G(2) = Y axis acceleration 
G(3) = Z axis acceleration 
l,G(i) supplied in tank triad 

0,G(i) supplied in inertial triad; Euler angles must 
be supplied to relate tank triad to this inertial 
triad . 

Euler angle permutation type if ITRIAD=0 


ITYPE 

= 



1 

2 

3 

4 5 

6 

7 

8 

9 

10 

11 

12 

axis 

of 

1st 

rotation 

X 

X 

X 

X Y 

Y 

Y 

Y 

z 

z 

z 

Z 

axis 

of 

2nd 

rotation 

Y 

Y 

z 

Z Z 

Z 

X 

X 

X 

X 

Y 

Y 

axis 

of 

3rd 

rotation 

Z 

X 

X 

Y X 

Y 

Y 

Z 

Y 

z 

z 

X 


ANG(l-3) = X,Y and Z Euler rotation angles as specified by 
ITYPE above (in degrees) 


EE?R0DUC''P,?1.ITY OF THE 
ORiCINA'' ■ '•'■.AIE POOR 
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presented near the end of this appendix. The primary data from 
program TANK used in program LAMPS3 consists of: 

1. elliptical surface coefficients (a,c): page 6 of the 

sample output; 

2. liquid mass: page 8 of the sample output; 

3. spherical coordinate location (R, 0, 0) of the initial 

liquid cm position: corrected values page 10 of the 

sample output. "Corrected" refers to the fact that 
the liquid cm has been adjusted to the"best fit" 
ellipsoidal constraint surface. 

Having run program TANK the input for LAMPS3 is simple to 
implement. Table A-2 delineates the coded input for the LAMPS3 
sample case. The user is referred to Reference 3, Volume II 
for a detailed definition of the input for DYNAMO. The input 
presented for LAMPS3 is identical in format to the DYNAMO in- 
put for nonlinear time response of interconnected rigid bodies. 
In Table A-2 the names of the input variables that change to 
reflect a specific tank/liquid system are shown in parenthesis; 
ie, (Rq)*, Table A-3 details the nature of these variables for 
the user. All other input is unchanging and reflects the 
mechanization of the three dimensional model. 

Sample output for LAMPS3 is presented at the end of this 
appendix. The user is directed to Reference 3 for a full 
definition of the output variables. At each integration time 
interval (as selected by the PR input parameter) the values of 
the state vector and its derivative, Y and YDT, are printed. 
There are 34 variables in each vector (as defined in equation 
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REPRODUCffilLnY OP THE 
ORIGINAL PAGE IS POOR 





DEN 410740 (6-74) 



FORMA 







FORMA 
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TABLE A- 3. INPUT VARIABLE DEFINITIONS FOR LAMPS3 


VARIABLE 


DESCRIPTION (FORMAT)^ 

RUNNO 

— 

Run number printed in page heading 1 3A6) 

User name printed in page heading J ’ ’ 

UNAME 


TITLE 1 

= 

Title card printed in page heading (12A6) 

TITLE2 

- 

Title card printed in page heading (12A6) 

COMMENT 

” 

Comment cards, last card = zeros in columns 
1 through 10 

9Y,0Z,9X 


Initial f3 *s for hinge Q2/P2 (Figure A- 1) 
0Y = 0 radians 
0Z = -(V2) + THETA radians 
0Y = 7T - PHI radians 

where: THETA = output of program TANK (if 



PHE=0° or PHL=180° use THETA= 



^/2 radians) 



PHI = output of program TANK 

^o 

= 

Start time , seconds 

A T 


Time interval for integration, seconds 

Te 

= 

End time for the program, seconds 

PR 


Print every PRth time interval 

PL 

= 

Plot every PLth time interval 

a ,c 
K 


Ellipsoidal surface coefficients from program 
TANK 

— 

Differentiation algorithm gain 

FX,FY,FZ 

= 

Forces applied to Body 1, see Figure A- 1 

MASFAC 

= 

Effective mass f ac tor 

Ro 

= 

Spherical coordinate distance to liquid cm at 

M 

= 

Liquid damping, based on 0 (used in equation 



(A- 11)) 

Ml 


Mass of the drop capsule (or rigid spacecraft) 
Body 1 

^XX’^yy 

= 

Principle moments of inertia of Body 1 

M3 

= 

Liquid mass (Body 3) from program TANK 

PTITLE 

= 

This portion of the plot titles may be changed 
to reflect problem being simulated 


^ Most of the input changes are to matrix inputs. These matrices 
are input via subroutine READ which is a FORMA subroutine. For 
details on subroutine READ see: R.L.Wohlen, "Synthesis of Dynamic 
Systems Using FORMA - Fortran ^trix Analysis". Martin Marietta 
Corp., MCR-71-75, Vol. IV, NAS8-25922, May 1971. Note: on the 
coding sheets (Table A-2) the allowable field for the numerical 
inputs is delineated by underlining, ie, ^ 0. (QZ)" , - 
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111-13). The specific variables in the state vector (Y) for 
our mechanization are as follows, (the variables in YDT are 
just the derivatives of the Y variables) : 


STATE VECTOR LOCATION VARIABLE 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 
, 33 

34 


z 

u 

V 

w 


“z 

u 

V 

w 


(Body 1 - Tank (Rigid Spacecraft) 


l-Body 2 - Dummy 


bBody 3 - Liquid 


/I's : 

I Body 1 to 
l-Inertial Triad 

I ^ * s : Body 2 to Body 1 

L j8 : Body 3 to Body 2 

iRehonomic Z 3/2 
ySpherical Coordinates 
of Liquid cm 

; j Liquid Forces 
; Exerted on the 
;L Tank (Body 1) 


The liquid forces exerted on the tank. Body 1, are printed in 
the YDT vector; locations 32, 33 and 34. 
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The spherical coordinates of the liquid cm position are 
also found in the YDT vector; locations 29 and 30 (9,0). The 
spherical coordinate, R, is found by: 

R = - Y(27) = Ro - Z 3/2 

As mentioned earlier, the success of this implementation was 
dependent on differentiating i . The error in this differen- 
tiation can be assessed by comparing Y(28),J'«^ , and YDT(27), 
^ 3 / 2 * ^oth of these values are equivalent to R. LAMPS3 plots 
these two variables for comparison. Along with the sample out- 
put presented at the end of this appendix are the plots gener- 
ated by LAMPS3. Plots of FY and FZ from the two dimensional 
model LAMPS are also presented for comparison with the liquid 
forces generated by LAMPS3. 

It will be noted that the two dimensional model, LAMPS, 
predicts slightly higher forces than LAMP S3 . This is a result 
of the differences in implementation of the kinematics between 
the two programs, LAMPS3 includes the force feedback onto 
Body 1 due to the moving fluid, which results in altered gross 
accelerations during the course of the simulation. Program 
LAMPS does not include this feedback. In LAMPS, accelerations 
are applied directly to the liquid mass as opposed to forces 
applied to Body 1 in LAMPS3 . The kinematics of LAMPS3 are 
correct for simulating actual spacecraf t/ liquid interactions. 
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SAMPLE OUTPUT PROGRAM TANK 



9IJN NO, TANKi:^ OATF 10FE76 “ARE NO, \ 

»I)N nv Rl OEPPV 

11,«*f,51 CLOCK TIME 

TEST I :? PHASE 1 5,?40 CPTtNE 

.10 "lOOrl to «EC. PPTTHE 
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SAMPLE OUTPUT PROGRAM LAMPS3 


.R!»K NO, TLM13 n«TE 0?FE76 PftGE NO. 1 

RUN BY' DYNAMO 

t?.O0.5O CLOCK TIME 

SIMULATION OF TLM TEST CASE 13 IPHASE II 3f.7«*0 SEC. CPTIMF 

transient LTQUin MOTION.. .PHIsO DEC, 25 7. FILL 29 SEC. PPTIHE 
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APPENDIX B - TEST RESULTS 


This appendix presents (Figures B-1 through B-12) the 
reduced test data for all 23 tests. Shown in these figures 
are plots of the tank Y and Z forces vs time. 

Table B-1 summarizes the test conditions for each case. 
Values of axial acceleration (AZI) were calculated based on 
test time, drop capsule and drag shield masses, and drop 
capsule travel distance. Lateral accelerations (AYI) were 
scaled from high speed photographs taken during testing. 
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TABLE B-1. AXIAL AND LATERAL ACCELERATION TABLE 


TEST # 

TANK 

CONFIGURATION 

ex 

(DEG) 

7o FILL 

TIME 

(SEC) 

AZI 

(IN/SEC^) 

AYI 

(IN/SEC^) 

1 

Baseline 

0.0 

25. 

0.0 

386.04 

0.0 





0.07 

-33.20 

0.0 





1.64 

-33.20 

0.0 

2 

Baseline 

45.0 

25. 

0.0 

386.04 

0.0 





0.08 

-32.25 

0.0 





1.66 

-32.25 

0.0 

3 

Baffled 

0.0 

10. 

0.0 

386.04 

9.72 





0.08 

-32.25 

9.72 





1.66 

-32.25 

9.72 

4 

Baf 

fled 

45.0 

10. 

0.0 

386.04 

10.12 






0.07 

-31.79 

10.12 






1.67 

-31.79 

10.12 

5 



90.0 

10. 

0.0 

386.04 

10.39 






0.09 

-30.89 

10.39 






1.69 

-30.89 

10.39 

6 



0.0 

25. 

0.0 

386.04 

10.43 






0.09 

-28.38 

10.43 






1.75 

-28.38 

10.43 

7 



45.0 

25. 

0.0 

386.04 

9.91 






0.07 

-28.78 

9.91 






1.74 

-28.78 

9.91 

8 



90.0 

25. 

0.0 

386.04 

7.78 






0.08 

-28.98 

7.78 






1.735 

-28.98 

7.78 

9 



0.0 

50. 

0.0 

386.04 

7.931 






0.08 

-27.68 

7.931 





1.768 

-27.68 

7.931 
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TABLE B-1. (continued) 


TANK 


Small 


TIME 

(SEC) 

AZI 

(IN/SEC^) 

AYI 

(IN/SEC^) 

0.0 

386.04 

6.53 

0.07 

-27.99 

6.53 

1.76 

-27.99 

6.53 

0.0 

386.04 

7.87 

0.07 

-26.12 

7.87 

1.81 

-26,12 

7.87 

0.0 

386.04 

8.73 

0.06 

-16.40 

8.73 

2.05 

-16.40 

8,73 

0.0 

386.04 

7.24 

0.08 

-16.21 

7.24 

2.06 

-16.21 

7.24 

0.0 

386.04 

7,77 

0.07 

-16.40 

7.77 

2.05 

-16.40 

7.77 

0.0 

386.04 

10.34 

0,07 

-28.78 

10.34 

1.74 

-28.78 

10.34 

0.0 

386.04 

8.990 

0.08 

-28.71 

8.990 

1.742 

-28.71 

8.990 

0.0 

386.04 

8.051 

0,07 

-28.585 

8.051 

1.745 

-28.585 

8.051 

0.0 

386.04 

8.63 

0.08 

-28.386 

8.63 

1.750 

-28.386 

8.63 



TABLE B“l. (continued) 


TIME 

(SEC) 

AZI „ 
(IN/SEC^) 

AYI 

(IN/SEC^) 

0.0 

386.04 

8.50 

0.06 

-16.498 

8.50 

2.045 

-16.498 

8.50 

0.0 

386.04 

8.52 

0.07 

-16.595 

8.52 

2.040 

-16.595 

8.52 

0.0 

386.04 

20.96 

0.08 

- 28.386 

20.96 

1.28 

- 28.386 

20.96 

0.0 

386.04 

21.12 

0.08 

- 28.386 

21.12 

1.295 

- 28.386 

21.12 

0.0 

386.04 

7.96 

0.08 

-27.606 

7.96 

1.77 

-27.606 

7.96 



